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1. INTRODUCTION 

The control of a bipedal robot is a challenging task due to the many degrees of freedom involved 
and the nonlinear and hard-to-stabilize dynamics [1], [2]. Bipedal robot structure is one of the most flexible 
forms of a walking robot. A bipedal robot has the same mechanism of action as a human and capable of 
walking in an environment containing uneven terrain, slope, stairs, obstacles, etc. [3]-[7]. 

Bipedal stability during walking motion is a critical fact in preventing the robot from falling down 
and causing the human or itself damages. The Zero Moment Point (ZMP) method is often used as a stability 
criterion for a walking biped robot [8]-[14]. We have found in the literature the use of force sensors placed 
under the feet. In [15] the authors proposed a Fuzzy control scheme in which the controlled variable is the 
coordinate compensation of the trunk only in the sagittal plane, used to move the measured ZMP to the 
desired ZMP, which is obtained from force sensors placed at the robot’s feet. 

Another related work [16] also proposes a cooperative control scheme to achieve lateral balance of a 
biped robot in the coronal plane, tested under disturbing conditions such as changes in the floor inclination 
and external forces applied to the body. The control system is composed of one master controller in charge of 
commanding the legs servomotors and two proportional Fuzzy controllers, one for each foot, which compute 
the Center of Pressure (CoP) under each foot and generate a control action to move the hips-ankles 
servomotors in order to return the CoP to a stable region. 

In [17] analyzed ZMP was analyzed from human walker wearing robot feet as shoes. Some 
researchers proposed a flexible shoe system for biped robots to an optimized energy consumption of the 
lateral plane motion [18]. Kitti Suwanratchatamanee et al.[19] proposed a haptic sensing foot system for 
humanoid robots. They investigate two different kinds of implementations: one is an active tactile sensing 
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technique to recognize a contacting ground slope. The other is to balance the robot body with one leg for 
human-robot interaction. 

Reference [20] proposed a fuzzy controller to realized biped walking on the slope with the desired 
ZMP trajectory. In [21] proposed a dynamic balance control (DBC), which includes Kalman filter (KF) and 
the fuzzy motion controller (FMC), was also designed to keep the body balance and make the biped walking 
following the desired ZMP reference. In addition, KF is utilized to estimate the system states and reduced the 
effect caused by noise. 

Stability is a vital criterion of bipedal robots in dynamic or static walking. The CoP located at the 
position under the foot bottom should affect bipedal walking stability. Thus the study of the CoP region for 
which a bipedal robot remains stable is a chalenge. Stable walking of a bipedal robot can be realized by 
balancing two types of force. First is the resultant force of gravity and the inertia produced by the motion of 
the robot. The other is the ground reaction force acting on both feet of the robot. In this paper, we developed 
stability walking bipedal robot in the second perspective. We refer the point associated with contact forces as 
CoP. 

In this paper, we propose a method to control bipedal robot to have a stable walking on uneven 
terrain. Through an intelligent control such as fuzzy logic control (FLC). The major contributions of this 
paper are follows: 1) development and improvement of walking stability of bipedal robot based on Force 
Sensitive Resistor (FSR) sensor; 2) development of a foot pad sensor by using FSR; 3) design a FLC for CoP 
region's stability and stable compliance control; and 4) presents the practical effectiveness of this control 
scheme by experiments with HuroEvolutionJR kid size bipedal robot. 

The remainder of this paper is organized as follows. In Section II, Strict definitions of CoP and ZMP 
are specified and body posture control bipedal robot also walking method of a bipedal robot. Section III 
elaborates the control system architectures especially on a fuzzy system for this work. Section IV presents 
several simulation results which show the effectiveness and merit of the proposed methods. Section V 
concludes the paper with remarks and suggestions for future works 


2. RESEARCH METHOD 
2.1. Center of pressure (CoP) and zero moment point (ZMP) definition 

The concept of ZMP [22] is introduced under the assumption that a robot walks on a flat floor and 
the frictional forces are strong enough to compensate the ground reaction forces tangential to the ground. 
ZMP, is the point on the surface of the foot where a resultant force R (Fr) can replace the force distribution 
shown in Figure 1. Mathematically, the ZMP can be calculated from a group of contact points P; for i=1,..., 
N with each force vector f; associated with the contact point. 


N 
Í: 
Dit : (1) 


In this definition, the ZMP can never leave the support polygon. If the floor is assumed horizontal, 
the torque reduces to at the ZMP. The position of CoP, which is a point on the ground plane, where the 
ground reaction force is applied, is often considered to be an equivalent of ZMP, since they coincide while 
the robot is balanced. Therefore, the actual position of ZMP can be found using force sensors located on the 
soles of a robot. With these sensors, the center of pressure can be calculated on the feet, and the ZMP can be 
directly measured. 


Foot Right Support Pollygon 


Fr 
Foot Left 
CoP 


Figure 1. Support polygon and CoP position based on the total sum of the contact forces FR acts 
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2.2. Determining CoP position 

We can calculate the resultant force by using the center of area (CoA) technique. In this paper, we 
developed four force sensitive resistor (FSR) sensors which were attached to the sole of the foot. There are 
two step calculations in order to know the CoP position which is related to the gait phases, i.e, single support 
phase (SSP) and double support phase (DSP). Figure 2 shows the CoP position when the bipedal robot on the 
SSP. The coordinates A, B, C, D are described by P4 (Xpa, Ypa), Ps (Xp, Yrs), Pc (Xpc, ypc) and Pp (Xpp, Yep). 
So, we can calculate the CoP position as follows for SSP phase as indicated in (2)-(3). 


_ faXpa + feXpe + foXpc + fpXpp 


fatte+feotfo (2) 


CoP 


_ faYpat fe ¥en t+ fodec t+ foypp 


a+ fet fe+tfo (3) 


YCoP 


Where fi, fe, fc, fp, are the force at A, B, C, D point achieved through FSR sensors and explained in (4)-(7). 


Xp, = FL—dxA; yp, = FW—dyA 


(4) 
Xpp = FL—dxB ; ypg = dyB (5) 
Xpc = AKC ; Ype = FW —dxC (6) 
Xpp =dxD ; Ypp = dxD (7) 


Foot Length ( FL ) 


Figure 2. CoP position when the bipedal robot on SSP phases 


Figure 3(a) shows DSP gait phases. The coordinate Al, B1, Cl, D1, A2, B2, C2 and D2 are 
described by PAI (xpa, Ypai), PB1 (Xppi, Yppi), PC1 (Xpci, Ypci), PD1 (Xpp1, Yppi), PA2 (Xpaz, Ypa2), PB2 
(Xpp2, Ypp2), PC2 (Xpc2, Ypc2) and PD2 (Xpp2, Ypp2). When the bipedal robot on stand still conditions, the ideal 
CoP position is depicted on Figure 3(b). In this condition the force sensing is distributed with same value at 
all point. The CoP position located in circle black. 
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Foot Left 


Foot Right 


(a) (b) 


Figure 3. (a) CoP position when the bipedal robot on DSP phases and (b) CoP position when the bipedal 
robot on stand still conditions 


2.3. Foot pad sensor development 

FSR is sensors that allow detecting physical pressure, squeezing and weight. The FSR is made of 
two layers separated by a spacer. The higher force will push the active element dots closer to the 
semiconductor and makes the resistance go down. In this paper, we used FSR from Interlink 402 model. The 
easiest way to measure FSR sensor is using a voltage divider technique as shown in Figure 5. The junction 
point resistor up with FSR resistor is connected to the analog input of a microcontroller such as an Arduino. 
Four FSR sensors are equipped at the corners of each foot as shown in Figure 4. In order to minimize impact 
and deformation, and also to distribute the repulsive force, the sole is composed of plate up, plate down and a 
rubber plate. The FSR sensors are fixed between the two plates 
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Figure 4. Foot pad sensor development architectures 


Suppose a FSR sensor has already set by a calibration processed, then it can proceed for further 
investigation. As mentioned in the datasheet, FSR data is characterized by a low frequency region. Therefore, 
our system compensates those sensor data by implementing low-pass filter (LPF). Theoretically, a LPF 
denote as in (8). 
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dv 
out 


dt (8) 


v.(t)-v (t) = RC 


out 


Let’s simplify the input and output at evenly spaced points in time, then we can derive simplicity as in (9). 


Y Ya A RC 
x =y, =RO— Tp t Ya Y: = aX, + (l-a)y,, 
A, RC+A, RC+A, (9) 
A, 1 . . . 
where @ , RC= , f. denotes as the cut-off frequency, A, is the sampling period and the 
RC+A, 2a f. 


smoothing factor œ lies in Osa@sl. 


rmm, 
memos 
» 


D 


Vout = 


Low Pass Filter Ri+R2 


Figure 5. The FSR sensors are fixed between the two plates and voltage divider configuration at FSR with 
low pass filtering 


2.4. Body posture control on bipedal robot 

There are two ways in order to control the body posture bipedal robot, i.e. inverse kinematics and 
inverse kinematic with adding some theta offset to each joint motor. In three dimensional motion, frontal 
plane (roll movement) and sagittal plane (pitch movement) as shown in Figure 6(a) are contained with ten 
joints servo motor each leg. Table 1 is described the frontal and sagittal plane joint names 

There are several parameters that will be used to adjust the posture of a bipedal robot as depicted in 
Figure 6(b). In this paper, we adjusting the hip_x_offset of our bipedal robot. The range operation of 
hip_x_offset is from -3 until 3 with cm units. 


Table 1. The Frontal and Sagitall plane joint names of Bipedal Robot 


Frontal plane (roll) Sagittal plane (pitch) 
OL left hip roll Ar left hip pitch 
O47 left ankle roll O37, left knee pitch 
OR right hip roll O51, left ankle pitch 
O4R right ankle roll OR right hip pitch 


OzR right knee pitch 
OsR right ankle pitch 
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Figure 6. (a) Sagittal and frontal plane in a bipedal robot (b) Hip_x_offset parameters 


2.5. Walking method of bipedal robot 

Based on the human way of walking, bipedal robot walks divide into two phases, namely single 
support phases (SSP) and double support phases (DSP) [23]. In the beginning of walking, both of the feet are 
in touch with the base, which called DSP as described in Figure 7. The overall locomotion control system can 
be designed as shown in Figure 8. At the beginning , the (XemdYemoOcma) Were commands that may be 
generated from an autonomous navigator or user's command, where (XgmaVendOcma) are relative translational 
landing position of a swing foot with respect to a support foot; Oena is the relative landing foot heading angle. 
Hence, the (XgnaVend9ena) Command is an omnidirectional walking command for bipedal robots that can be 
used to process footprint placement trajectory planning [24]. 


Walking Direction 


Single support 


Double support Single support Single support 


Figure 7. Human walking gait phase 
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Figure 8. Bipedal robot locomotion control system 


To achieve the trajectory of planes locomotion, the end position of each foot is also required. The 
trajectory of the swing foot is generated from a cycloid curve which can be indicated in (10)-(12). It is noted 
that (x, Ys Zs) are the position of the swing foot; Length, Height and Shift are the desired strike length, the 
strike height and the shift distance; p is the time percentage of the period time when the foot reaches the 
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highest position; T,is the period time. To get the real robot rotation based on the desired rotation, (13) is 
utilized to do so, where Ty, is the rotation time constant. As a consequence, the gait trajectory can be obtained 
as the swing foot trajectory and the CoM trajectory. The inverse kinematics (IK) can be applied to the gait 
trajectory to find the joint angles of the bipedal leg model. The joint commands were sent to the servomotors. 


Length t t 
(at) = 2 in(2 ; O<t<sT 
x, ( ) ae [27 T. sin( TT J] k (10) 
ee ee OSS 
y(t) = m sin(27 ; Stsi, 
s m T T : (11) 
Height t . t 
[27 sin(2z )] O<t< pT, 
2a pT, pT, 
Z(t) = s z m 
Height t— t— 
Heien- LEM (22 Pe ae E pT, <t <T, 
27 a- p)T, a- p)T, i i (12) 
aoe 
; 2 
hiP rotation Taa o a 
T,Ad-2Toc) (13) 


2.6. Fuzzy logic controller (FLC) 

Originally advocated by Zadeh [25] and Mamdani and Assilian [26], fuzzy has become a mean of 
collecting human knowledge and experience and dealing with uncertainties in the control system process as 
shown in Figure 9(a). The design of fuzzy logic controller generally has the following steps: 

a. Design the membership function for fuzzification input and output variables. 

b. Implement the fuzzy inference by a series of IF- THEN rules. 

c. Inference engine derives a conclusion from the facts and rules contained in the knowledge base using 
various human expert techniques. 


d. Design of defuzification for converting fuzzy data into numerical data that can be transmitted to the 
control equipment. 
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fa 
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Figure 9. (a) Fuzzy architectures and (b) Fuzzy inference system (FIS) from MATLAB 


MATLAB® Fuzzy Logic Toolbox is used to aid in FLC design as shown in Figure 9(b). The 
toolbox contains functions, graphical user interfaces and data structures that allow the user to quickly design, 
test, simulate and modify a fuzzy inference system. In this paper using Mamdani fuzzy rules as inference 
rules. This method is often also known as Max-Min method. This method was introduced by Ebrahim 
Mamdani in 1975. 
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To obtain the output fuzzy, it takes four stages: 

1. Fuzzification is the process of transforming crisp values into degree of membership for linguistic terms of 
fuzzy sets. In this paper using three types of analysis of fuzzy set membership in fuzzification process, 
there are listed as follows. 

a. The upgraded fuzzy membership function is described in Figure 10(a) and Equation (14). 


0 x<a 
=4 (x-a) a<x<b 
(b = a) (14) 
b. The downward fuzzy membership function is described in Figure 10(b) and Equation (15). 
0 x<a 
f= )(x-a) a<x<b 
b-a) (15) 


c. The membership function of the triangular fuzzy is described in Figure 10(c) and equation (16). 


0-7) a<x<b 
(ba) 
p= 
(x-a) 
= b<x<c 


(16) 


(a) (b) (c) 


Figure 10. Type of membership function in Fuzzy; (a) upgraded fuzzy, (b) downward fuzzy, 
(c) triangular fuzzy 


2. The fuzzy implication is to evaluate the consequent part of each rule. Among various implication methods 
available in the literature, Mamdani implication method is selected. After the inputs have been fuzzified, 
and fuzzy logic control knows the degree to which part of the antecedent rule is satisfied. The equation of 
Mamdani method to verify the fuzzy implication is in (17). 


Hy (k) = max [ min { u, (k), HM, (error(i),derror( D] an 


3. Generating the rules for a fuzzy control system is often the most difficult step in the design process. It 
usually requires some expert knowledge of the plant dynamics. There are three methods used in fuzzy 
system inference, namely: max, additive and probabilistic OR. In the max method, the fuzzy set solution 
is obtained by retrieving the rule maximum value, then using it to modify the fuzzy region, and applying 
it to the output using the OR (union) operator. 

4. We proposed weight average (WA) for defuzzification stage [27]. The WA calculation is described in 
(18) using three singletons membership function 


WA = u(k1)xkl+u(k2)xk2+ u(k3)xk3 
u(k1) + u(k2) + u(k3) (18) 
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3. FUZZY SYSTEM ARCHITECTURES 
3.1. Fuzzy logic controller on CoP region stability 

The FLC is one of the effective methods to control a system which is hard to model clearly, like a 
bipedal robot system. The control algorithm can be interpreted and formulated by expert knowledge. Before 
designing the FLC, the desired control goal should be defined first. In this section, we described the FLC on 
CoP region’s stability. The basic idea is to keep balance the CoP position of a bipedal robot when walks. In 
SSP, the CoP position should be located inside the sole of the foot. On the other hand, it should be located 
between the two feet in the DSP as illustrated in Figure 11. 


Rule Base 


Y 
” is ig - B 
j ae, Mechanism Aes 
and Fuzzification Inference —m Defuzzification | _“ _» 
o! 
y-axis 
x 
(a) (b) 


Figure 11. (a) CoP position of bipedal robot when walks at DSP phase and (b) the FLC design in stability 
region 
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Figure 12. The x-axis and y-axis of margin stability region 
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Figure 13. Fuzzification (a) x-axis, (b) y-axis membership function of CoP region's stability 


From the result of determining CoP position, we are able to define the area of CoP position as the 
input of our FLC system. In x-axis and y-axis of margin stability region as shown in Figure 12, we can define 
three memberships of fuzzy sets, N (Negative), Z (Zero) and P (Positive). After a fuzzification process of a 
fuzzy set x-axis and y-axis, we constructed rule based FLC in order to keep the CoP position still inside of 
stability regions. The relationship between CoP position on the x-axis and y-axis will be implicated by using 
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Table 2. The fuzzification x-axis and y-axis are using triangle membership function with 3 fuzzy sets as 
shown in Figure 13. The length value of membership function for x-axis and y-axis are 12 cm and 20 cm as 
indicated the total foot length of the robot. 


NS ZE PS 


p 


3 o 3 (cm) 


Defuzzification with 3 set singleton 


Figure 14. defuzzification membership function of CoP region's stability 


Table 2. Rule base for CoP region's Stability 


x-axis CoP Position 


freien) ai | rs | 
0 Z 0 


c 
S 
> Negative (N) 
£ 
3 
S Zero (2) 
š N Zz P 
~ 
Positive (P) 0 z 0 


The next stage is defuzzification. The membership function of this stage is described in Figure 14. 
The result of defuzzification on this case is the length of hip_x_offset. This length will use for adjusting the 
robot body posture tends to lean backward or forwards. When the robot walking on uneven terrain as shown 
in Figure 15, the feedback of CoP position will be evaluated by FLC, and then output this controller will 


adjust the hip_x_offset 
a «// 


hip_x_offset 


| 
After CoP Before CoP 
feedback feedback 


Figure 15. Body posture changes after CoP feedback 


3.2. Fuzzy logic on stable compliance control 

In this section, we explained the compliance control with implementing using FLC. The basic idea 
is to keep balance the bipedal robot when any force applied on the feet. For an example, when we pushed the 
robot from the back side. The measurement of force sensor on front side feet more than back as illustrated in 
Figure 16. 
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Figure 16. Motion after and before pushing from back direction 


In compliance control, we also measure and consider the CoP position of bipedal robot. In the ideal 
condition, when the bipedal robot stand-up on DSP, the position of CoP will in the center (black circle) as 
shown in Figure 17(a). The stability margin, red line rectangle is used to avoid and boundaries the CoP 
position which still inside of stability margin. By knowing the real position of CoP, we will know the real 
condition of bipedal robot. For example, when the CoP position at condition 1 (see Figure 17(b), we can 
measure the average force on points B and D at foot left. If the force point B and D at foot left more than our 
force limitation which is mean more force on left side foot. Thus, we can define the error as in (19). 


error = positive pitchbody limitation— pitchbody robot measurement (19) 


Furthermore, if we know the error, we could calculate the delta error as indicated in (20). 


deltaerror| = error —lasterror 
(20) 


Foot Length (FL) 


Foot Width 
(Fw) 


Foot Right 


(a) (b) 
Figure 17. (a) CoP position of bipedal robot when stand up at DSP phase (b) four categories of CoP position 


and its condition 


After computing the absolute error and the absolute delta error as indicated in Equations (19) and 
(20), these data are taken as input variables of the designed FLC for condition 1 case as shown in Figure 18. 
The ouput process of defuzzifiction is offset joint (@). 
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Figure 18. The FLC design in stable compliance control 


The specification of fuzzy design are explained as follows: 1) fuzzification error and delta error using 
triangle membership function with 5 fuzzy sets as shown in Figure 19, 2) rule base evaluation using MacVicar 
Whelan style inference as shown in Table 3, and 3) defuzzification using 5 set singletons as shown in 
Figure 20. 


Fuzzification membership function error Fuzzification membership function derror 


Figure 19. Fuzzification membership function on compliance control 
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Figure 20. Defuzzification membership function on compliance control 


Table 3. Rule base for Compliance Control 


error 
Very Small Large Very Large 
(NB) (Ps) (PB) 
Very Small 
(NB) NB NB NS NS z 
Small 
d (NS) NB NS NS Z PS 
error ea TNS NS z PS PS 
Large 
(Ps) NS Z PS PS PB 
oe EZ PS PS PB PB 


(PB) 


After the defuzzification process, we get one output (Ø) that will be used as value offset joint on 
ankle, knee, and hip as explained in Figure 21(a) and Figure 21(b). 
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Figure 21. (a) Motion before and after at compliance control for condition and (b) Motion before and after at 
compliance control for condition 2 


Figure 21(a) is described when the bipedal robot in condition 1. In this situation, the value output 
from defuzzification is theta offset ¢ that will be used to adjust the joint offset of roll for hip-left, ankle-left 


and also right. The formulation of this case is indicated in (21). 


hip _roll _left = -(@,, + ¢) 
ankle _ roll _left =—(@,, + ¢) 
hip _roll _ right = -(0,, + $) 
ankle _ roll _right = -(0,, + 9) 


For condition 3, the formula is indicated in (22). 


hip _ roll _left = (0,, +9) 
ankle _roll _left = (0,, + 9) 
hip _roll _right = (0,, +) 
ankle _ roll _ right =(0,, + 9) 


(21) 


(22) 


Figure 21(b) is described the bipedal robot in condition 2. In this situation, the value output from 
defuzzification is theta offset ¢ that will be used to adjust the joint offset of pitch for the hip-left, knee-left, 


ankle-left and also right. The formulation of this case is indicated in (23). 


hip _ pitch _left =k, *(6, + g) 
hip _ pitch _right = k, *(0, + g) 
knee _ pitch _left =k,*(0,, -¢) 
knee _ pitch _right =k, * (O — $) 


ankle _ pitch _left =k, *(0,, -¢) 
ankle _ pitch _ right = k,*(6,, -¢) 


For condition 4, the formula is indicated in (24). 


(23) 
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hip _ pitch _left =—(k,*(0,, + ¢)) 
hip _ pitch _ right = -(k, #( 6. + ¢)) 
knee _ pitch _ left = -(k, *(8,- $)) 
knee _ pitch _ right = -(k, * (O = ¢)) 
ankle _ pitch _left =-(k,*(0,, -¢)) 


ankle _ pitch _right =—(k, *(0,, -)) (24) 


4. RESULTS AND ANALYSIS 
4.1. Bipedal platform 

In this paper, we used HuroEvolutionJR kid-size bipedal robot from Medical Robot Laboratory 
(MRL) National Taiwan University of Science and Technology. The HuroEvolutionJR is designed as an 
eighteen degree-of-freedom (DOF); where 12 DOF joints are desired for two legs, 4 DOF joints are desired 
for two arms, and 2 DOF joints are desired for the head. The height and the total weight including the 
batteries are 380 mm and 2900 gram, respectively. The mechanical structure of is shown in Figure 22(a). 

The electronic design of HuroEvolutionJR is shown in Figure 22(b), we used Arduino DUE ARM 
(Advanced RISC Machine) 32-bit microcontroller, namely Atmel SAM3X8E ARM Cortex-M3 CPU with 84 
MHz speed of frequency clocks. The Arduino DUE is used for motion controller with deals for trajectory 
patterns such as Linear Inverted Pendulum Model (LIPM), Inertial Measurement Unit (IMU) sensor 
processing, and FSR sensor processing. This motion controller will have communication data with personal 
computers (PC) and remote control device using serial data communication. 


Arduino 8 bit 


Host PC 
GUI - Processing 


(a) (b) 


Figure 22. (a) Physical of the HuroEvolution’® (b) Electronic architectures 


4.2. Simulation result 

In order to verify the CoP equations on SSP and DSP, simulation by using MATLAB was applied 
before implement it on the bipedal robot. Two condition simulations (SSP or DSP) is shown in Figure 23. In 
this simulation, the force value is begins from 0 (indicated that no pressure) until 255 (indicated that high 
pressure). Futhermore, we also developed graphical user interface (GUI) to diplay foot pad sensor condition 
in realtime to evaluate the CoP conditions during static and dynamic walking as shown in Figure 24. This 
GUI is developed by Processing IDE. 
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FSR calculation during SSP FSR calculation during DSP 
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Figure 23. CoP position during SSP with fA=192, fB=128, fC=50, fD=90 [left] and CoP position during DSP 
with displacement=3 cm, fA1=125, fB1=56, f{C1=238, fD1=43, fA2=0, fB2=0, fC2=12, fD2=80 [right] 
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Figure 24. CoP position when simulate in GUI, CoP position changes to right 


4.3. Experimental result 

When the bipedal robot moves, the CoP trajectory was used to evaluate the CoP position which still 
inside of the CoP margin stability range or inside of support polygon as shown in Figure 26(a). For analyzing 
the result, a single walking cycle was divided into four parts: In Figure 25 shown that phase A-B is the first 
lifting of the right leg. Phase B-C is SSP on the right foot. Phase C-D is SSP on the left foot. And D-E phase 
is SSP on the right foot, this process will be repeated. 

There are three cases in order to verify the compliance control based on CoP condition and the force 
when applied in feet of the robot. Case one, when the CoP position changing to the left side of the robot as 
shown in Figure 28. Case two, when the CoP position changing to the front side of the robot as depicted in 
Figure 29(a) to Figure 29(b). Case three, when the CoP position changing to the back side of the robot as 
shown in Figure 29(c) to Figure 29(d). The fuzzy logic that applied in these cases is able to achieve for 
keeping balance the body posture of the robot. The response of error, delta error and the output of 
defuzzification will be shown in Figure 26(b). The parameter of membership function fuzzification error and 
delta error are tuned by experiment with GUI shown in Figure 27. Finally, Figure 30 is shown the conditions 
of a bipedal robot when walks in uneven terrain. All video experiments were upload at YouTube 
channel in [28]. 
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DSP SSP right foot DSP SSP left foot 
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Figure 25. The video capture during the robot’s walking 


CoP Trajectory y-axis Stable Compliance Control 
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Figure 26. (a) The signals response error, derror and (b) Ouput fuzzy when CoP position changing to left side 
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Figure 27. Fuzzy logic tools GUI to evaluate the output of fuzzy 
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Motion before pushing Motion after pushing 
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Figure 28. Stable compliance control when CoP position changing to left side; (a) ilustration (b) real 
condition record 


Motion before pushing Motion after pushing 


© 
Motion after pushing Motion before pushing 


Figure 29. Stable compliance control when CoP position changing to front side; (a) ilustration 
(b) real condition record and Stable compliance control when CoP position changing to back side; 
(c) ilustration (d) real condition record 
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Figure 30. Stable compliance control when CoP position changing to front side 


5. CONCLUSION 

In this paper, a method of controlling bipedal robots to walk stably is proposed based on FSR 
sensor. Fuzzy logic control (FLC) was used for two cases. The first was for stable compliance control and the 
second was CoP region's stability. In compliance control mode, the designed FLC accepts inputs from error 
and delta error of force limitation and generate theta offset that used as feedback to adjust the body posture of 
a bipedal robot. In CoP region's stability mode, the designed fuzzy controller accepts inputs from x-axis and 
y-axis of CoP position and changes the robot body posture tends to lean backwards or forwards. 

Eventually, several experiments showed that the designed of FLC able to make the robot walks 
stable on uneven terrain. For future works, this system can be added with another sensor like vision. With 
vision perception, the robot can recognize the terrain condition in environments so when the robot walks in 
an uneven terrain, the recognized terrain condition would provide useful information. 
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